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We consider the possibility that a large fraction of the ultrahigh energy cosmic rays are decay 
products of Z bosons which were produced in the scattering of ultrahigh energy cosmic neutrinos on 
cosmological relic neutrinos. We compare the observed ultrahigh energy cosmic ray spectrum with 
the one predicted in the above Z-burst scenario and determine the required mass of the heaviest relic 
neutrino as well as the necessary ultrahigh energy cosmic neutrino flux via a maximum likelihood 
analysis. We show that the value of the neutrino mass obtained in this way is fairly robust against 
variations in presently unknown quantities, like the amount of neutrino clustering, the universal 
radio background, and the extragalactic magnetic field, within their anticipated uncertainties. Much 
stronger systematics arises from different possible assumptions about the diffuse background of 
ordinary cosmic rays from unresolved astrophysical sources. In the most plausible case that these 
ordinary cosmic rays are protons of extragalactic origin, one is lead to a required neutrino mass in the 
range 0.08 eV < m„ < 1.3 eV at the 68 % confidence level. This range narrows down considerably 
if a particular universal radio background is assumed, e.g. to 0.08 eV < m„ < 0.40 eV for a large 
one. The required flux of ultrahigh energy cosmic neutrinos near the resonant energy should be 
detected in the near future by AMANDA, RICE, and the Pierre Auger Observatory, otherwise the 
Z-burst scenario will be ruled out. 

PACS numbers: 14.60.Pq, 98.70.Sa, 95.85.Ry, 95.35.+d 



I. INTRODUCTION 

Big bang cosmology predicts the existence of a back- 
ground gas of free photons and neutrinos. The measured 
cosmic microwave background (CMB) radiation supports 
the applicability of standard cosmology back to photon 
decoupling which occured approximately one hundred 
thousand years after the big bang. The relic neutrinos, 
on the other hand, have decoupled when the universe had 
a temperature of one MeV and an age of just one second. 
Thus, a measurement of the relic neutrinos, with a pre- 
dicted average number density of 



56 cm" 3 . 



(1) 



CMB 



per light (m Vi <C 1 MeV) neutrino species i, would pro- 
vide a new window to the early universe. Their pre- 
dicted number density is comparable to the one of the mi- 
crowave photons. However, since neutrinos interact only 
weakly, the relic neutrinos have not yet been detected 
directly in laboratory experiments p], |^, ||, ||, [|, |6[ Q . 

Recently, an indirect detection possibility for relic neu- 
trinos has been discussed || ||]. It is based on so-called 
Z-bursts resulting from the resonant annihilation of ul- 
trahigh energy cosmic neutrinos (UHEC^s) with relic 
neutrinos into Z bosons Jk], |l[ || (cf. Fig. |l|). On 



resonance, the corresponding cross section is enhanced 
by several orders of magnitudes. If neutrinos have non- 
vanishing masses m Vi - for which there is rather con- 
vincing evidence in view of the apparent observation of 
neutrino oscillations |l3[ |l4i |15l [lq, |j7|, |l8|, |lS|, |20|, 21, 

H, n n || n, nrHlllMiplI U m SJ B 

- the respective resonance energies, in the rest system of 
the relic neutrinos, correspond to 



ET = 



Ml 
2 m Vi 



= 4.2 • 10 21 eV 



1 eV 



(2) 



with Mz denoting the mass of the Z boson. These 



uhec v 



o 



10 71 

2 nucleons 
n Tt " + — - e 



20 y 



,v,v 



* fodor@pms2.elte.hu 



5audor.katz@dcsy.de 



$ andrcas.ringwakl'idcsy.dej 



FIG. 1: Illustration of a Z-burst resulting from the resonant 
annihilation of an ultrahigh energy cosmic neutrino on a relic 
(anti-)neutrino. 
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resonance energies are, for neutrino masses of 0(1) 
eV, remarkably close to the energies of the highest en- 
ergy cosmic rays observed at Earth by collaborations 
such as AGASA §7|l,_Fly's Eye || || [|J, Haverah 

(for a review, 




that the ultra- 



Park f4lj,l2j, HiReS |3|, and Yakutsk 
see Ref. ]45|]). Indeed, it was argued 
high energy cosmic rays (UHECRs) above the predicted 
Grciscn-Zatsepin-Kuzmin (GZK) cutoff 46, [47) around 
4 • 10 19 eV are mainly protons (and, maybe, photons) 
from Z decay. In this way, one possibly also solves one of 
the outstanding problems of ultrahigh energy cosmic ray 
physics fl48| , namely the apparent observation of cosmic 
rays with energies above the GZK cutoff, in an elegant 
and economical way without invoking new physics be- 
yond the Standard Model, except for neutrino masses. 
The GZK puzzle hinges on the fact that nucleons with 
super-GZK energies have a short attenuation length of 
about 50 Mpc, due to inelastic interactions with the cos- 
mic microwave background, while plausible astrophysi- 
cal sources for those energetic particles are much farther 
away [|9| Ultrahigh energy neutrinos produced at 
cosmological distances, on the other hand, can reach the 
GZK zone unattenuated and their resonant annihilation 
on the relic neutrinos could just result in the observed 
cosmic rays of the highest energies. Moreover, at present 
this annihilation process may be the only way to detect 
the relic neutrino background, a basic ingredient of our 
cosmological picture. 

The Z-burst hypothesis for the ultrahigh energy cosmic 

fs was discussed in many p ap e rs fill, p2] |53| , |5~| |55] |57] . 
|58| , |59| , |60| , pi] , |62] , p3| , |6J, |65j |66|. The required 
UHEC^ fluxes were estimated in Ref. |n] for different 
spectral indices. In Ref. J52|, particle spectra were deter- 
mined numerically for case studies which supported the 
Z-burst scenario. The effect of possible lepton asymme- 
tries was studied in Ref. [Q . In Ref. , the analysis of 
the Z-burst mechanism was advocated as one of the few 
possibilities for an absolute neutrino mass determination. 

In the present paper, we present the details of (and 
extend) our recent quantitative investigation of the Z- 
burst scenario pi] (see also Refs. [62, pij), where we 



have determined the required mass of the heaviest relic 
neutrino as well as the necessary ultrahigh energy cosmic 
neutrino flux via a maximum likelihood analysis. But 
before we start this enterprise, we shall briefly review 
the current knowledge of the absolute scale of neutrino 
masses. 

Neutrinos almost certainly have non-vanishing masses 
and mix. This follows from the apparent observation 
of neutrino oscillations whose evidence is compelling for 
atmospheric neutrinos jl 3| , fl 4[ Jl 5j , ilj 17], strong for 
solar neutrinos ||, [Tgfcj! ^lplglT^rEl , and so- 
far unconfirmed for neutrinos produced in the labora- 
tory and studied by the LSND collaboration and oth- 
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How- 



ever, neutrino oscillations are sensitive only to the mass 



(squared) splittings Am 
vidual masses, (m Vi >) 



i?,., not to the indi- 



Only a lower bound on the mass of the heaviest neutrino 
can be derived from these observations, e.g. 



m V3 > jAml tm > 0.04 eV 



(3) 



from the atmospheric mass splitting in a three neutrino 
flavour scenario and 



> 



A? 



LSND 



> 0.4 eV 



(4) 



from the still allowed value of the LSND mass splitting 
in a four flavour scenario, respectively (for a recent com- 
pilation of global analyses on mass splitting and mixing 
parameters, see for example Refs. |i7] ^8] |6j| and refer- 
ences therein). 

For an investigation of the absolute scale of neutrino 
masses one has to exploit different types of experiments, 
such as the search for mass imprints in the endpoint spec- 
trum of tritium beta ((3) decay or the search for neu- 
trinoless double beta (0v2(3) decay. At present, these 
direct kincmatical measurements of neutrino masses pro- 
vide only upper limits, e.g. 



m/3 



/^|E/ ei | 2 m2. < 2.2-=-2.5eV (95%C.L.), (5) 



with U being the leptonic mixing matrix, for the effective 
mass measured in tritium (3 decay [(70|, fn], [72], (73|, and 



< 0.33 -4-1.35 eV (90%C.L.)(6) 



for the Majorana neutrino mass parameter |74|, [7g, 
[77| , [78[ appearing in 0v2(3 decay. The recently reported 
evidence for neutrinoless double beta decay and corre- 
spondingly deduced parameter range Q 



(m„) = (0.11 -4- 0.56) eV (95 % C.L.) 



(7) 



has been seriously challenged by Refs. (gCj, |8lj. Com- 
bining the experimental constraints from oscillations and 
from tritium (3 decay, one infers upper bounds on the 
mass of the heaviest neutrino |B2l, 



m U3 < x jra\ + Am^ < 2.5 eV . 



in a three flavour, and 



in. 



< 



m \ + Ato lsnd ^ 3 - 



eV. 



(8) 



(9) 



in a four flavour scenario, respectively. For further re- 
cent investigations of the relationship between oscillation 
phenomena, (3 decay and Qv2f3 decay, in particular with 
respect to the neutrino mass spectrum, one may also con- 
sult Refs. fl83|, |84|, I 



36 87 



V3 > m„ 2 > m Ul , themselves. 



and references 

cited therein. 

Further information on the absolute scale of neutrino 
masses can be obtained through cosmological and astro- 
physical considerations. Neutrinos in the 0.1-^1 eV mass 



3 



range have cosmological implications since they represent 
a non-negligible part of dark matter. This gives the op- 
portunity to put upper limits on neutrino masses from 
cosmology Analyses of galaxy clustering, in- 

cluding recent CMB measurements and other cosmolog- 
ical constraints, give an upper bound 



< 1.8 -7- 4.4 eV 



(10) 



on the sum of the neutrino masses ||| ^4 , 95 , 9(| . From 
the spread of arrival times of neutrinos from supernova 
SN 1987A, coupled with the measured neutrino energies, 
a time-of-flight limit of mp < 23 eV can be derived |97], 
|98|, which, however, is not competitive with the direct 
limit (||). 

It is extremely welcome that the Z-burst scenario opens 
a further and timely window to the absolute neutrino 
mass scale, since the opportunities to determine this cru- 
cial quantity are rare at present. In addition, its verifica- 
tion would give us an indirect detection of the so elusive 
relic neutrinos and, finally, offer an explanation of the 
origin of the highest energy cosmic rays. 

The organization of the present paper is as follows. In 
Section O we describe our determination of the proton 
and photon spectra at Earth, which originate from Z- 
bursts taking place mainly at extragalactic distances. We 
discuss in detail the main ingredients in the predictions 
of the spectra, such as the details of Z production and 
hadronic decay, the propagation of nucleons and photons 
through the diffuse extragalactic photon background, the 
diffuse flux of ultrahigh energy cosmic neutrinos, and the 
relic neutrino number density, as well as their anticipated 
uncertainties. The comparison of the Z-burst spectra 
with the observed ultrahigh energy cosmic ray spectrum 



is presented in Section III. The required absolute neu- 



trino mass, as well as the necessary ultrahigh energy cos- 
mic neutrino flux, are determined via a maximum likeli- 
hood analysis for various assumptions about the nature of 
the background of ordinary cosmic rays from unresolved 
astrophysical sources and variations of the diffuse extra- 
galactic photon background notably in the radio band. 
On the basis of these studies we find that the Z-burst 
determinations of the mass of the heaviest neutrino as 
well as of the neutrino flux are fairly robust. We find a 
required neutrino mass range of 0.08 eV < m v < 1.3 eV 
at the 68 % confidence level, if the background of ordi- 
nary cosmic rays is of extragalactic origin. This range 
narrows down considerably if a particular universal radio 
background is assumed, e.g. to 0.08 eV < m v < 0.40 eV 



for a large one. In Section IV we discuss the implications 



of our findings for future laboratory studies such as the 
tritium beta decay and neutrinoless double beta decay, 
as well as for astrophysical and cosmological neutrino in- 
vestigations. This section contains also our conclusions. 



II. Z-BURST SPECTRA 

Our comparison of the Z-burst scenario with the ob- 
served UHECR spectrum proceeds as follows. First, we 
determine the probability of Z production as a function 
of the distance from Earth. Secondly, we exploit col- 
lider experiments to derive the energy distribution of the 
produced protons and photons in the laboratory (lab) 
system. Thirdly, we consider the propagation of the pro- 
tons and photons, i. e. we determine their energy losses 
due to pion and/or e + e~ production through scattering 
on the diffuse extragalactic background photons and due 
to their redshift. The last step is the comparison of the 
predicted and the observed spectra and the extraction of 
the required mass of the heaviest relic neutrino and of 
the necessary UHECV flux. 

Our prediction of the differential proton flux, i.e. the 
number of protons arriving at Earth with energy E per 
units of energy (E), of area (A), of time (t) and of solid 
angle (17), 



Fp\z(E) 



dN„ 



p\Z 



dE dA dt dn ' 

from Z-bursts can be summarized as 



F, 



p\Z 



CO JXma 

(E)=Y,J dE P J 



dr 



(11) 



(12) 



dE Vi F n {E Vi , r) n^(r) 



■DC 

J dE 9i F Pi (E^.,r)n Vi (r 



Xfr/jPi (s) Br(Z — * hadrons) 



dN p+n 
dE n 



x(-)—P p (r,E p ;E), 

with the following important building blocks: the 
UHECzy fluxes F Vi [E Vi , r) at the energies E Vi (» E T v f) 
and at the distance r of Z production to Earth, the 
number density n Vi (r) of the relic neutrinos, the Z pro- 
duction cross section <r Vt i> t (s) at centre-of-mass (cm) en- 
ergy squared s = 2 m Vi E Vi , the branching ratio Br(Z — > 
hadrons), the energy distribution dN p+n / dE p of the pro- 
duced protons (and neutrons) with energy E p , and the 
probability P p (r, E p ; E) that a proton created at a dis- 
tance r with energy E p arrives at Earth above the thresh- 
old energy E. 

A similar expression as Eq. ([l2]) holds for the differen- 
tial photon flux from Z-bursts, 



F-y\z(E) = E/ dS 7 / 



dr 



(13) 
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OO 

J dE Vi F Vi {E v .,r)n 0i {r) 



+ / AE 9i F^{E^,r)n Vz {r) 
o 

dA 

xcr^.p. (s) Br(Z — > hadrons) 7 

dE~, 



Here, the photon propagation function P 7 (r, E 7 ] E) gives 
the expected number of detected photons above the 
threshold energy E if one photon started from a distance 
of r with energy E 1 . Note, that the P 1 function has a 
different interpretation than the P- function of protons. 
This arises from the fact that the number of photons - in 
distinction to the number of protons - is not conserved 
during their propagation. 

The building blocks related to Z-production and decay 
~~ GviVi i the hadronic branching ratio, and the momentum 
distributions dNi / dEi - are very well, and the propaga- 
tion functions Pi are fairly well determined, whereas the 
first two ingredients, the flux of UHECVs, F Vi (E Vi ,r), 
and the radial distribution of the relic neutrino number 
density n Ui (r), are much less accurately known. In the 
following we shall discuss all these ingredients in detail. 



A. Z production and decay 

At LEP and SLC millions of Z bosons were produced 
and their decays analyzed with extreme high accuracy. 
Due to the large statistics, the uncertainties of our anal- 
ysis related to Z decay turned out to be negligible. 

Let us start with a review of the Standard Model 
neutrino annihilation cross section. The s-channel Z- 
exchange annihilation cross section into any fermion anti- 
fermion (//) pair is given by (see e.g. Refs. @ ||) 



G 



47T (s - 



Ml? 



M%Y% 



N eS (s) , (14) 



where s is the cm energy squared, T z is the total width 
of the Z boson, and N e g is the effective number of anni- 
hilation channels, 



iV eff (s)=^0(s-4m))x 
/ 

2 

X - rif (l - 8 T 3 f qf sin 2 9 W + 8 qj sin 4 



(15) 



Here the sum is over all fermions with m/ < \/s/2, with 
charge qf (in units of the proton charge), isospin T3 f 
(1/2 for u, c and neutrinos; —1/2 for d, s, b and negatively 
charged leptons), and rif = 1(3) for leptons (quarks (<?)). 



With sin 2 6 W = 0.23147(16) for the sin 2 of the effective 
Weinberg angle and G F = 1.16639(1) x 1(T 5 GeV~ 2 for 
the Fermi coupling constant ]100| , formula ( H ) gives, at 
the Z-mass, 

a(uii/i -> Z* -> all qq) \ S=M 2 = 314.9 nb , (16) 

aiyiPi -» Z* -> all //) | s=M s = 455.6 nb , (17) 

with a branching ratio aiyiVi — * Z* — > all qq)j a(yiVi — > 
Z* —> all //) | s= m| = 0.6912, in good agreement with 
the experimental result flOOfl , 

(18) 



Br(Z -> hadrons) = (69.89 ± 0.07) % . 



Later we shall exploit the following simplification 
which arises due to the fact that the cross section ( |l4| ) 
is sharply peaked at the resonance cm energy squared 
s = M§ . Correspondingly, it acts essentially like a 5- 
function in the integration over the energies E Vi of the 
incident neutrinos in Eqs. ([l2|) and (|l^), and we can as- 
sume that the UHECV fluxes are constant in the relevant 
energy region. Thus, introducing the energy-averaged 
annihilation cross section [l(| , 

Knn) = J = 2 ^ ^ G F = 40.4 nb , (19) 

which is the effective cross section for all neutrinos within 
1/2 6E Z ™/E Z ™ = Tz/Mz = 2.7% of their peak annihila- 
tion energy, we can write 

OO 

J dE Vi F Vi {E Vi ) a ViVi (s = 2 m Vi K, ) (20) 
o 

~ e™ F Vi (K7) W- 

In view of the expected rapid decrease of th e UHECz; 
flux at increasing energies (cf. Section |ll C| ), we ne- 
glect t-channel W- and Z-exchange annihilation pro- 
cesses. On resonance, the s-channel Z-exchange processes 
completely overwhelm them. Fairly above the resonant 
energies, on the other hand, they eventually dominate 
but are probably unobservable due to the lack of an ap- 
preciable UHECV flux. 

Next, we turn to the energy distribution of the protons 
and photons in Z decay. Wc combined existing published 
and some improved unpublished data on the momentum 
distribution, 



_ P P 



Pbc 



(21) 



of protons (p) (plu s antiprotons (p)) in Z decays [101 



102] , |0§ [Uj3 |l05| , see Fig. | (top). The experimental 
data, ranging down to x w 8 ■ 10~ 3 , were combined with 
the predictions from the modified leading logarithmic ap- 
proximation (MLLA) 



(TM [107], |10§ at low x. The p + p 
multiplicity is (N p ) = J dxV p {x) = 1.04 ± 0.04 in the 
hadronic channel 11100(1 . 
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FIG. 2: Momentum distributions in hadronic 



Combined data from Refs. 101 



102, 10E 



104 



Z de cays. Top: 
on proton 



105 



(plus antiproton) momentum distribution (solid) , normalized 
to (N p ) — 1.04, and charged pion momentum distributions 
(dotted), normalized to (N v ±) = 16.99. Bottom: Distribution 
of protons ("p"; solid), photons ("7"; dotted) and electrons 
("e"; dashed) in the lab system, in which the target relic 
neutrino is at rest. 



In the cm system of the Z production the angular 
distribution of the hadrons is determined by the spin 
1/2 of the primary quarks and thus proportional to 
1+w 2 = 1+cos 2 9 (here 9 is the angle betwee n the incom- 
ing neutrinos and the outgoing hadrons (cf. |109|] )). The 
energy distribution of the produced protons with energy 
E p entering the Z-burst spectrum (|l2|), 



dE p ~ E v dy ~ E v Jp[V) ' 



(22) 



with y = 2E p /E u , is finally obtained after a Lorentz 
transformation from the cm system to the lab system, 



Qp(v) = ^lj^ dw ( 1+w2 ) 



(23) 



1 



±y - w^Jy 2 - (1 - w 2 )(2m p /M z y 



^y 2 ~{\-w 2 )(2m p /M z Y 



1 — w 2 

f-wy ± y/y 2 - (1 - w 2 ){2m p /M z ) 2 



p 



1 



where m p is the proton mass. The first line comes from 
the angular distribution, the second line is the Jacobian 
and the third one is the momentum distribution at the 
inverted momentum. The scaled energy distribution Q p , 
as defined in Eq. ( p2| ) and given by Eq. (23), is displayed 
in Fig. U (bottom). 

Neutrons produced in Z decays will decay and end up 
as UHECR protons. They are taken into account accord- 
ing to 



Qp+n{y) = 1 



(Nn) 



Qpiv) 



(24) 



where the neutron (n) (plus antineutron (n)) multiplic- 
ity, ( N n ) = J dy Q n (y) , is ~ 4% smaller than the pro- 
ton's fuct 

Photons are produced in hadronic Z decays via frag- 
mentation into neutral pions, Z^tt° + X^2^ + X (cf. 
Fig. [I]) . The corresponding scaled energy distribution in 
the lab system, defined analogously to Eq. (E3) , reads 



dw ■ 



2.y 



w 



(25) 



where the scaled energy distribution Q 7r o(y), with y — 
2E n o/E„, is given by Eq. fl23|), with m p — > m^a and 
V p — > T^o, the momentum distribution of pions in 
hadronic Z decay. For the latter distribution, we take 
the me asured one of charged p ions V^± from hadronic Z 
decay |o| |o| [l0| |l0|, |o| (cf. Fig. | (top)), normal- 
ized such that (N y ) = J^dyQ^iy) = 20.97 p0| . Note, 
that the distributions Qi, i = p, 7, presented in Fig. || 
(bottom), compare favorably with the ones prese nted in 
Ref. |36| based on the event generator PYTHIA [111]. 

Electrons (and positrons) from hadronic Z decay are 
also relevant for the development of electromagnetic cas- 
cades. They stem from decays of secondary charged pi- 
ons, Z—>n ± +X—>e ± +X (cf. Fig. |l|), and their scaled 
energy distribution in terms of y = 2E e /E v reads 



-1 p2 

2 e ±(y) = I dw I dx 



1 



1 jo xw + y/ x 2 + (2 m e /m w ) 2 

2y 



(26) 



xw + \/ x 2 + (2r7i e /m 7r ) 2 



Ve(x) . 



where V e is the momentum distribution of the electrons 
in the rest system of the charged pion. The energy dis- 
tribution Q e ± is also displayed in Fig. |^ (bottom). 
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B. Propagation of nucleons and photons through 
the diffuse extragalactic photon background 

The cosmic microwave background is known to a high 
accuracy. It plays the key role in the determination of 
the probability P p (r, E p ; E) that a proton created at a 
distance r with energy E v arrives at Earth above the 



threshold energy E, suggested in Ref. [ 1 1 2 1 and deter- 
mined for a wide range of parameters in Ref. [ 1 13 1 . The 
propagation function P p takes into account the fact that 
protons of extragalactic (EG) origin and energies above 
w 4 ■ 10 19 eV lose a large fraction of their energy due 
to pion and e + e~ production through scattering on the 
CMB and due to their redshift Q fuj. In the present 
analysis we have included new results for P p (r, E p ; E) 
which include now variations in the cosmological param- 
eters, in extension to the already published form |113| 
exploited in Ref. foil . In our analysis we go, according 
to 



dz = -(l + z)H(z)dr/c, 



(27) 



out to distances i?max (cf. (|12|)) corresponding to redshift 
Zmax = 2 (cf. Ref. [114 1). We use the expression 



h 2 {z) = hI [n M (i + z) 3 + n A ] 



(28) 



for the relation of the Hubble expansion rate at redshift 
z to the present one. Uncertainties of the latte r, Hn = h 



100 km/s/Mpc, with h = (0.71 ± 0.07) x J 
included. In Eq. (28), S1m and fi A , with Q,m 



1.15 
95 



[ lOOp , are 

-n A = i, 

are the present matter and vacuum energy densities in 
terms of the critical density. As default values we choose 
Qm = 0.3 and fi A = 0.7, as favored today. Our results 
turn out to be pretty insensitive to the precise values of 
the cosmological parameters. 

The simulations needed for the computation of 
P p (r, E p ] E) require large computer power. We have used 
a farm of personal computers (PCs) consisting of 128 par- 
allel 1.7 GHz Pentium 4 pro cesso rs - normally devoted 
to QCD lattice calculations [ 1 1 5 1 - and have exploited 
about 3 hours CPU time to determine P p (r, E p ; E) in a 
range of r < 4000 Mpc, 10 18 eV < E p < 10 26 eV, and 
10 18 eV < E < 10 26 eV, for each fixed value of the cos- 
mological parameters. The simulation was carried out 
in small (10 kpc) steps in r. For each step, the sta- 
tistical energy losses due to pion/e + e~ production and 



redshift are taken into account 1 1 13 ] . In this connection, 
the advantage of our formulation of the Z-burst spec- 
trum (cf. (|12|)) in terms of the probability P p (r, E p ; E) 
becomes evident. We have to determine the latter only 
once and for all. Without the use of P p (r, E p ; E), we 
would have to perform a simulation for any variation of 
the input spectrum, notably for any change in the neu- 
trino mass. Our maximum likelihood analysis, involving 
the neutrino mass as a free parameter, would thus require 
excessive computer power - on the order of 300 hours on 
the above PC farm for fixed cosmological parameters. 
Since P p (r, E p ; E) is of universal usage, we have decided 
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FIG. 3: Top: The intensity spectrum of the diffuse extra- 
galactic photon background at redshift z = from Ref. [116] 
(solid). Different estimates from Ref. [118] of the universal 



radio background (URB) are also indicated: high URB (long- 
dashed) and moderate URB (short-dashed). Bottom: Pho- 
ton energy attenuation length la at z — c orresponding to 
the photon background shown above Jll6| (solid). Varia- 
tions of In aris ing from different assumptions about the URB 
(cf. Ref. 1 118] and top) are also indicated: high URB (long- 
dashed) and moderate URB (short-dashed). Shown is fur- 
thermore the energy attenuation length for electrons due to 
synchrotron radiation (long-dashed-short-dashed), for differ- 
ent magnitudes of the extragalactic magnetic fields. 



to make the corresponding numerical data for the proba- 
bility distribution (— )dP p (r, E p ; E)/dE available for the 
public via the World- Wide- Web URL 



http :/ / www, desy.de/ ~uhecr 



The determination of the photon propagation function 
Pj(r, Ey, E), entering the photon flux prediction (p"3|), 
was done as follows. In distinction to the case of the pro- 
ton propagation function, we used here the continuous 
energy loss (CEL) approximation which largely simplifies 
the work and reduces the required computer resources - 
we have estimated the necessary CPU time for a full sim- 
ulation on the PC farm mentioned above to 1000 hours 
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per fixed set of the cosmological parameters. In the CEL 
approximation, the energy (and number) of the detected 
photons is a unique function of the initial energy and dis- 
tance, and statistical fluctuations are neglected. A full 
simulation of the photon propagation function will be the 
subject of a later work. 

The processes that are taken into account are pair pro- 
duction on the diffuse extragalactic photon background 
(cf. Fig. | (top)), double pair production and inverse 
Compton scattering of the produced pairs. We comment 
also on synchrotron radiation in a possible extragalac- 
tic magnetic field (EGMF). For the energy attenuation 
length of the photons due t o the se processes, we exp loited 
the values quoted in Ref. |116|] (see also Ref. [117Q and 
the further ones presented in Fig. [| (bottom) which incor- 
porate various assumptions abo ut th e poorely known uni- 
versal radio (URB) (from Ref. and infrared (IRB) 
backgrounds. We shall analyse later the dependence of 
the neutrino mass and other fit parameters on these vari- 
ations. Note , that, in view of the recent URB estimates 
in Ref. |118[ , the ones presented in Fig. || (top)), which 
are based on Ref. [US], can be referred to as "minimal" 
URB. 

The computation of the photon propagation function 
P 7 (r, E-y ; E) was carried out in the following way. The 
energy attenuation of photons in the CEL approximation 
was calculated according to 



dE = -E 



dr 



dz 



h{E) 1 + 2 



(29) 



where l z {E) = (1 + z)~ 3 l (E(l + z)) is the energy atten- 
uation length at redshift z. The number of photons was 
assumed to be constant at ultrahigh energies > 10 18 eV, 
due to the small inelasticities in this energy range. Below, 
it was increased in a way to maintain energy conservation 
(except for the redshift contribution): 



-iV-y 



l z {E) 



(30) 



The Pj{r, E~ l \ E) function was then obtained by integra- 
tion of these equations. In the ultrahigh energy region - 
which is most relevant for us since we perform our fit to 
the cosmic ray data there - the approximation described 
above gives the photon flux quite reliable, while at lower 
energies it yields an upper bound. 



C. UHECi/ fluxes 

Presently unknown ingredients in the evaluation of 
the Z-burst spectra (|l2] ) and (|l3|) are the differential 
fluxes F v . of ul t rahi gh energy cosmic neutrinos (see e. g. 
Refs. [120, 121, 122 1 for recent reviews). Present exper- 
imental upper li mits on th ese fluxes are rather poor (cf. 
Fig. | and Refs. f|§ p9| ). 

What are the theoretical expectations for diffuse 
UHEO fluxes? More or less guaranteed are the so- 
called cosmogenic neutrinos which are produced when 
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FIG. 4: Upper limits on differential neutrino fluxes in the ul- 
trahigh energy regime. Sho wn ar e experimental upper limits 
on F Ve + F Pe from Fly's Eye |l23) and on J] Q = e ( F » a + ) 
from the Gold stone lunar ultrahigh energy neutrino expe- 
riment GLUE [124 1, as well as th eore ti cal upper limits on 



F v „ + Fp fro m "visible" ("WB" 



125 



126 [) and "hidden" 



("MPR" [127 1) hadronic astrophysical sources. 



ultrahigh energy cosmic protons scatter inelastically off 
the cosmic microwave background radiation ^7| in 
processes such as — > A — > nn + , where the produced 
pions subsequently decay |130|, 131 L The s e flu xes (for 



U7j , [132] |133| ) repre 



recent estimates, see Refs. [[12| , 
sent reasonable lower limits, but turn out to be insuf- 
ficient for the Z-burst scenario. Recently, theoretical 
upper limits on the ultrahigh energy c osmi c neutrino 



27). Per con- 
visible" hadronic as- 



flux have been given in Refs. [125, 126 
struction, the upper limit from 
trophysical sources, i. e. from those sources which are 
transparent to ultrahigh energy cosmic protons and neu- 
trons, is of the order of the cosmogenic neutrino flux 
and shown in Fig. | ("WB"; cf. Refs. (l25[ jl2l| ). Also 
shown in this figure ( "MPR" ) is the much larger upper 
limit from "hidden" hadronic astrophysical sources, i. e. 
from those sources from wh ich only photons and neu- 
trinos can escape [127, 134 1. Even larger fluxes at ul- 



trahigh energies may arise if the hadronic astrophysical 
sources emit photons only in the sub-MeV region - thus 
evading the "MPR" bound in Fig. ^ - or if the neutri- 
nos are produced via the decay of superheavy relic par- 
" |13§ |13§ [137], [L3§ p|, pTC " 



tides 



for which 



also the fragmentation function of the decay is of ma- 
jor interest |l4l| , pT2] |l43] , |144[ pT5| , pTr| , or topological 
defects @ p7| . 

In this situation of insufficient knowledge, we take the 
following approach concerning the flux of ultrahigh en- 
ergy cosmic neutrinos, F Vi (E Vi , r). It is assumed to have 
the form 



F Vi {E Vi , r) = F Vi (E Vi , 0) (1 + z) a , 



(31) 



where z is the redshift and where a characterizes the cos- 
mological source evolution (see also Refs. |l^, 52, p5[). 
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The flux at zero redshift, F Vi {E u „Q) = F Vi {E Vi ), is left 
open. For hadronic astrophysical sources it is expected 
to fall off power-like, F Vi {E v ^) oc E~p , 7>1, at high 
energies. Due to this fact and because of the strong reso- 
nance peaks in the ViVi annihilation cross section (IT3) at 
the resonance energies (||), the Z- burst rate will be only 
sensitive to the flux at the resonant energy of the heaviest 
neutrino. Of course, the latter may be nearly degenerate 
with the other neutrino mass eigenstates, m Vi w m„, as 
it is the case for m„ 3 > 0.1 eV in a three flavour scenario. 
Correspondingly, our later fit to the UHECR data will be 
sensitive only to 



F* 



E 



[F Vi (E™) + F^E™ 



(32) 



where the sum extends over the number of mass eigen- 
states which are quasi-degenerate with the heaviest neu- 
trino. Note, finally, that, independently of the produc- 
tion mechanism, neutrino oscillations result in a uniform 
F Vi mixture for the different mass eigenstates i. 



D. Neutrino number density 

The dependence of the relic neutrino number density 
n Vi on the distance r is treated in the following way. 

The question is whether there is remarkable clustering 
of the relic neutrinos within the local GZK zone of about 
50 Mpc. It is known that the density distribution of relic 
neutrinos as hot dark matter follows the tot al m ass distri- 

1491. In fact, 



bution; however, with less clustering [148, 
for m„ < 1 eV, one expects pretty much that the neu trin o 
number density equals the big bang prediction (Q) Jl50[ . 
To take above facts into account, the shape of the n Vi (r) 
distribution is varied, for distances below 100 Mpc, be- 
tween the standard cosmological homogeneous case (|l|) 
and that of the total mass dis tribution obtained from pe- 
culiar velocity measurements [ 151 , 152f (cf. Fig. H (top)). 



These peculiar measurements suggest relative overdensi- 
ties of at most a factor f v = 2-^3, depending on the 
grid spacing (cf. Fig. || (bottom)). A relative overden- 
sity f v = 10 2 4- 10 4 in our neighbourhood, as it was as- 
sumed in earlier investigations of the Z-burst hypothe- 
sis H |[ |5l], ^2], |5^, m, seems unlikely in view of these 
data. Our quantitative results turned out to be rather 
insensitive to the variations of the overdensities within 
the considered range. 

For scales larger than 100 Mpc the relic neutrino den- 
sity is taken according to the big bang cosmology predic- 
tion, n Vi = 56 • (1 + z) 3 cm -3 . Possible uniform neutrino 
density enhancements due to eventual lepton asymme- 
tries, as advocated in Ref. |54|], are negligible in view 
of the rece nt, very s tringent bounds on the neutrino de- 



generacies [153, 154 1 . In any case, such uniform enhance- 



ments change only the normalization of the Z-burst spec- 
tra and have no effect on their shape. Therefore, they 
could possibly milder the required UHECi' flux, but the 
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FIG. 5: Top: Mass density fluctuation field 5 along the 
supergalact ic p lane as obtained from peculiar velocity mea- 
surements 1 151]. Shown are contours in intervals of S = 0.2, 
surface maps on a grid of spacing 500 kms -1 , correspond- 
ing to 5/i -1 Mpc, with the height proportional to S, and 
contrast maps. One recognizes some well-known structures 
in the nearby volume such as the Great Attractor at super- 
galactic coordinates (SGX ~ -2000 kms" 1 , SGY ~ -500 
kms -1 ), the Perseus-Pisces complex (SGX ~ 6000 kms , 

SGY < 1000 kms -1 ), and the large void (SGX ~ 2500 

kms -1 , SGY ~ kms -1 ) in between. Bottom: Mass density 
fluctuation field obtained from above data, averaged over all 
directions, for h = 0.71. The overdensities at around 20 and 
80 Mpc reflect the Great Attractor and the Perseus-Pisces 
complex, respectively. 



value of the neutrino mass inferred from the Z-burst sce- 
nario is not affected by them. 
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III. DETERMINATION OF THE REQUIRED 
NEUTRINO MASS AND THE NECESSARY 
UHECi/ FLUX 

A. Generalities 

The predicted spectra of protons and photons from Z- 
bursts, ([l^) and (fi"3|), can now be compared with the 
observed UHECR spectrum (cf. Fig. ||). Our analysis 
includes published UHECR data of AGASA @, Fly's 
Haverah Park gl], ft|, and HiReS g|, 



Eye 1 38 



as well as unpublished one from the World Wide Web 
pages of the experiments on 17/03/01 (for a review, see 
Ref. (4^]). Due to normalization difficulties we did not 
use the Yakutsk |44j results. We shall take into account 
the fact that above 4 • 10 19 eV less than 50 % of the 
cosmic rays can be photons _at the 95 % confidence level 
(C.L.) ■ 



156 



155| (see also Refs. 
As usual, each logarithmic unit between log(i?/eV) = 
18 and log(i?/eV) = 26 is divided into ten bins. The 
predicted number of UHECR events in a bin is taken as 



N(i) 



where £ 



fEi+i 

£ \ dE [F plbkd (E) 

JEi 



F p(+J)IZ (E)] , (33) 
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Hi 



■ s • sr is the total exposure (es- 



timated from the highest energy events and the corre- 
sponding fluxes) and where Ei — 10( 18+J / 10 ) eV is the 
lower bound of the i th energy bin. The first term in 



Eq. ( 



F, 



p|bkd 



corresponds to the diffuse background 



of ordinary cosmic rays from unresolved astrophysical 
sources. Below the GZK cutoff, it should have the usual 



and experimentally observed power-law form 45 1. The 
second term represents the sum of the proton and photon 
spectra, F p \ z + F^\ z , Eqs. ( |l2| ) and (|l3|), from Z-bursts. 

The separation of the flux into two terms (one from 
the power-law background and one from the Z-burst) is 
physically well motivated. The power-law part below the 
GZK cutoff is confirmed experimentally and, for extra- 
galactic sources, it should suffer from the GZK effect. In 
the Z-burst scenario, cosmic rays are coming from an- 
other independent source (Z-bursts), too. What we ob- 
serve is the sum of the two. As the detailed fits in the next 
section will show, the flux from Z-bursts is much smaller 
in the low energy region than the flux of the power-law 
background. Correspondingly, the low energy part of the 
spectrum (between 10 18 5 and 10 193 eV) has very little 
influence on the Z-burst fit parameters, notably on the 
neutrino mass. 

We shall study several possibilities for the background 
term F p | bkd . 

The first one is based on the assumption that the dif- 
fuse background of ordinary cosmic rays even above the 
GZK cutoff at 4 • 10 19 eV consists of protons which are 
produced in our neighborhood - within our galactic halo 
or at least within the GZK zone of about 50 Mpc. We 
shall refer to this possible cosmic ray background in the 
following as the "halo background". We should note, 
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FIG. 6: The available UHECR data with their error 
bars and the best fits (long-dashed) from ordinary cosmic 
ray protons originating from our local neighbourhood ("halo 
background"; top) and originating from diffuse extragalactic 
sources ( "EG background" ; bottom) , respectively. For the lat- 
ter case, the bump at 4 • 10 19 eV represents protons injected 
at high energies and accumulated just above the GZK cutoff 
due to their energy losses. The predicted fall-off for energies 
above 4 ■ 10 19 eV can be observed. 



however, that there are no apparent suitable astrophysi- 
cal sources within the GZK distance to which the highest 
energy events point (4^, - a fact which is difficult to 
reconcile with the halo background model. 

For the halo background, no GZK attenuation is in- 
cluded. The spectrum is assumed to have the usual 
power-law behavior which describes the data well for 
smaller energies J|7], |]| (cf. Fig. ^ (top)), 



£ F, 



7 P \bkd(E; A, (3) = 

a / e 

leV V leV / 



(34) 



(Halo bk'd) . 



A similar background model, however with a fixed power- 
law index (3 = 3.23, as found by AGASA to fit the data 
from 10 17 - 6 ^ 10 19 eV (cf. Table V in Ref.jp), has been 
exploited in a recent Z-burst simulation |66]| T 

In the second case, we assume that the diffuse back- 
ground of ordinary cosmic rays comes from protons 
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which originate from uniformly distributed, extragalactic 
sources - hence we name it "extragalactic background" . 
In view of the observed distribution of arrival directions 
of UHECRs, this assumption seems to be phenomeno- 
logically more realistic than the halo background model. 
Though there are some peculiarly clustered events, the 
overall distribut ion at present statistics seems to be prac- 



tically uniform [157, 15S, 159 



The extragalactic background suffers of course from 
GZK attenuation. Correspondingly, we take the above 
power-law ([54]), A ■ E~@, as an injection spectrum and 
take its modification due to the interactions with the 
CMB photons into account with the help of the proton 
propagation function P p (r, E p ;E), 

£F plbkd (E;A,(3) = 



dE p I dr (l + z(r)Y 
o 



A ( E, 



-13 



leV V leV , 
dP p (r,E p ;E) 



(35) 



dE 



(EG bk'd) 



The predicted spectrum of the extragactic background 
protons shows an accumulation at around the GZK scale 
4- 10 19 eV and a sharp drop beyond (see Fig. || (bottom)). 

As far as the Z-burst spectra entering our predic- 
tion |3^) are concerned, we proceed as follows. We shall 
mainly concentrate on the case where there is only one 
relevant neutrino mass scale m v , either because there 
are three neutrino types with nearly degenerate neutrino 
masses, m Vi ~ ?n u , or there is one neutrino which is 
much heavier than the other ones such that the contri- 
bution of the latter to the cosmic ray spectrum can be 
neglected since the corresponding resonance energies are 
much larger and the UHECV fluxes are expected to fall 
with increasing energy. Therefore, we fit only one neu- 
trino mass parameter m v . In this case, the spectra jl2] ) 
and ( |l3| ) for protons and photons from Z-bursts can be 
written as 



£F i[z (E;B,m u ) 



D 



4m„ 



Qi(y = 



dr (1 



z(r)) 3+a S n (r) 



Am u Ei 



(36) 



x(-) 



dPi(r,Ej;E) 
dE 



P,7; 



where a is the cosmological evolution parameter, 5 n (r) 
is the mass density fluctuation field (cf. Fig. g (bot- 
tom)), normalized to one, and Qi are the boosted mo- 
mentum distributions from hadronic Z decay, normal- 
ized to (N p+n ) = 2.04, for i — p, and to (iV 7 ) = 



2 (Ar 7r o) + (7V 7r ±) = 37, for i = 7. We are left here with two 
fit parameters, the mass m„ of the heaviest neutrino and 
the overall normalization B, which may be expressed, on 
account of Eqs. ( pfj| ) and (32), in terms of the original 
quantities entering Eqs. (12) and (O), as 



§ = B r(Z 



hadrons) i? max (n„ t ) (ia„„) E™ s F* es .(37) 



Note, that the neglection of finite width effects, 
1/2 5E™/ E T u f = Tz/Mz = 2.7%, in our implementa- 
tion (^6|) of the Z-burst spectra is perfectly adequate in 
view of the relative errors > 50 % which we will find later 
from our fits. 

The expectation value for the number of events in a bin 
is given by Eq. ([33]). To determine the most probable 
value for m v . we use the maximum likelihood method 



and minimize [ 160 1 the x 2 {PiA,B,m Vj ), 



X 2 = (38) 
lo s(w)= 260 

2 [N(i] - N Q (i) + N a (i) In (N (i)/N(i))] , 

1 °g(w)= 1 8-5 

where N Q (i) is the total number of observed events in 
the i th bin. Since the Z-burst scenario results in a quite 
small flux for lower energies, we take the lower bound 
just below the "ankle": E min = 10 185 eV. Our re- 
sults are insensitive to the definition of the upper end 
(the flux is extremely small there) for which we choose 
log(-E max /eV) = 26. The uncertainties of the measured 
energies are about 30% which is one bin. By means of 
a Monte Carlo analysis, we take these uncertainties into 
account and include the corresponding variations in our 
final error estimates. 

For comparison with recent work on the Z-burst sce- 
nario [Ea, we have done also fits with no background 



component, -F p |bkd = in Eq. (33), and a larger lower 



end, log(i? m in/eV) = 19.4 20.0. Such a scenario, with 
E m i n Eqzk ~ 4- 10 19 eV, was advocated in Ref. |65) as 
appropriate to explain all UHECR data above the GZK 
cutoff by the Z-burst model and to attribute all UHECR 
data below the cutoff to ordinary astrophysical extra- 
galactic sources. 



B. Fit results 

Qualitatively, our analysis can be understood as fol- 
lows. In the Z-burst scenario a small relic neutrino mass 
needs a large incident neutrino energy _E™ S (|2|) in order to 
produce a Z. Large E]f s results in a large Lorentz boost, 
thus large E p resp. E 1 . In this way the shape of the de- 
tected energy (E) spectrum determines the mass of the 
relic neutrino. The sum of the necessary UHECi^ fluxes 
FJ CS , on the other hand, is determined by the over-all 
normalization B. 



11 



TABLE I: Results of fits, for a strong UHE7 attenuation 
(h = 0.71, tl M = 0.3, fi A = 0.7,z max = 2). Top; Assuming 
a halo UHECR background according to Eq. (B4j). Middle- 
Assuming an extragalactic UHECR background according to 
Eq. (fj5|). Bottom: Assuming no UHECR background above 
.Emm, for different values of the lower end _B m in of the fit 
(a = 0). 



Halo UHECR background + strong UHE7 attenuation 


a 


m„ [eV] 


Xmin 


A 


B 





-3 

3 


A zl«+ 2 ' 22 ( 4 - 80 ) 
1. 64(2. 88) 
& c+l. 73(4.03) 
°- W -1.34(2.32) 
n r, +1.45(3.30) 
Z - J± -1.05(1.80) 


15.90 
15.64 
15.53 


1.46 ■ 10 43 
1.62 ■ 10 43 
1.65 ■ 10 43 


1049 
770 
551 


3.110 
3.111 
3.111 


EG UHECR background + strong UHE7 attenuation 


a 


m„ [eV] 


Xmin 


A 


B 





-3 

3 


n 9f) +0.20(0.63) 
U - ZU -0.11(0.18) 
n 9f) +0.19(0.61) 
U -^ U -0.12(0.18) 
„ „ n +0.19(0.59) 
U - ZU -0.11(0.17) 


25.82 
26.41 
26.89 


5.00 ■ 10 31 
5.98 ■ 10 31 
7.23 ■ 10 31 


150 
144 
142 


2.465 
2.466 
2.467 


No UHECR background + strong UHE7 attenuation 


log^min/eV) 


m„ [eV] 


Xmin 


A 


B 


P 


19.4 
19.5 
19.6 
19.7 
19.8 
19.9 
20.0 


9 9 o+0.64(1.46) 
Z ' Z -0.58(1.06) 
1 01+0.63(1.44) 
i " 3i -0.53(0.80) 
n pr+0.67(1.62) 
U - OJ -0.31(0.55) 

n zLn + - 32 (°- 87 ) 

U-^-0. 16(0.27) 
n -9+0.41(1.25) 
u -^ z _0.18(0.29) 
n 7fi +1.06(2.50) 
u '' "-0.39(0.58) 
7? +1.49(3.47) 
' -1.01(1.47) 


21.81 
16.01 
14.80 
8.03 
7.99 
5.52 
2.68 




1251 
846 
670 
445 
460 
733 
2021 





Our fitting procedure involves four parameters: /?, A, B 
and m v . The minimum of the x 2 (A A, B, m„) function 
is Xmin a t TO f min which is the most probable value for 
the mass, whereas the 1 a (68%) confidence interval for 
77i„ is determined by 



X 2 (/3', A', B', m v ) = Xo(m») = x mi „ + 1 



(39) 



Here f3' , A', B' are defined in such a way that the \ 2 
function is minimized in /3, A and B, at fixed m v . 

As already mentioned, presently there is no evidence 
that the observed highest energy cosmic rays are photons. 
Let us start therefore with the assumption (cf . Ref . (6lJ ) 
that the ultrahigh energy photons from Z-bursts can be 
neglected in the fit in comparison to the protons. This is 
certainly true for a sufficiently large universal radio back- 
ground, ejj. on the level of the maximal one estimated 
in Ref. [11? ] and/or for a sufficiently strong extragalactic 
magnetic field O(10~ 9 ) G. We shall refer to this scenario 
in the following as "strong" UHE7 attenuation. Our best 
fits to the observed data from this scenario (cf. Table |) 
can be seen in Fig. for evolution parameter a = 0. 

We find a neutrino mass of m v = 3.46^'g^2 '32) e ^ ^ or 
the case that the UHECR background protons are of 
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FIG. 7: The available UHECR data with their error bars 
and the best fits from Z-bursts, for a strong UHE7 at- 
tenuation such that the Z-burst photons can be neglected 
(a = 0,h = 0.71, Q M = 0.3, fi A = 0.7, ;w = 2). Top: 
Best fit for the case of a halo background (solid line). The 
bump around 4 • 10 19 eV is mainly due to the Z-burst pro- 
tons (dash-dotted), whereas the almost horizontal contribu- 
tion (long-dashed) is the first, power-law-like term of Eq. ([33]). 
Middle: The case of an "extragalactic" UHECR background. 
The first bump at 4 • 10 19 eV represents protons produced at 
high energies and accumulated just above the GZK cutoff due 
to their energy losses. The bump at 2 • 10 21 eV is a remnant 
of the Z-burst energy. The long-dashed line shows the con- 
tribution of the power-law-like spectrum with the GZK effect 
included. Bottom: The case of no UHECR background above 
log(£ min /eV) = 19.7. 



eV, if they are of 



12 



n 40 +a32 (°- 87 ) eV 

U - W -0. 16(0.27) eV ' 



extragalactic type fl35|), and m 

there are are no background protons above 10 19 7 eV, re- 
spectively. The first numbers are the 1 a, the numbers 
in the brackets are the 2 a errors. This gives an abso- 
lute lower bound on the mass of the heaviest neutrino of 
m„ > 0.02 eV at the 95% C.L., which is comparable to 
the one obtained from the atmospheric mass splitting in 
a three flavour scenario, Eq. (||). 

The surprisingly small uncert ainties are based on the 



if ^ 



X analysis described in Section [II A. The inclusion of 



the already mentioned 30 % uncertainties in the observed 
energies by a Monte Carlo analysis increases the error 
bars by about 10%. Note, that the relative errors in 
the extragalactic and in the no background cases are of 
the same order. This shows that the smallness of these 
errors does not originate from the low energy part of the 
background component. 

The fits are rather good: for 21 non- vanishing bins and 
4 fitted parameters they can be as low as Xmin = 15-64 
and Xmin = 25.82 in the halo and the extralactic back- 
ground case, respectively, whereas in the no background 
case, for E mm — 10 19 7 eV, we have 9 bins with 2 fit- 
ted parameters and a x min = 8.03, see Table |[ In the 
latter case, however, which was advocated strongly in 
Ref. [p5| , a remarkable dependence of the fitted mass on 
the value of E m [ n is observed. The ± 2 a fits are shown 
in Fig. ||. As it should be, the spread is small in the re- 
gion where there are data, whereas it can be quite large 
in the presently unexplored ultrahigh energy regime. Fi- 
nally, let us mention that, both in the "halo" as well as 
in the extragalactic background cases, the x 2 fits with- 
out a Z-burst component (cf. Fig. |^) are far worse: we 
find Xmin = 41.37 in the former, and Xmin = 38.25 in 
the latter, for 21 non-vanishing bins and 2 fitted param- 
eters. This finding suggests that the power-law back- 
ground terms alone cannot describe the data. 

In addition to the case of strong UHE7 attenuation, 
corresponding to a large universal radio background or 
a large extragalactic magnetic field, let us consider now 
the case of a "minimal" URB. Here, we a ssume a van- 
ishing EGMF and exploit, as in Ref. [116], a universal 



radio backgr ound on the level of the one from Ref. 1 1 19 
(cf. Section |f| and Fig. " 



(top)). The corresponding 
fits, for our three background scenarios, can be seen in 
Fig. § and Table || Note, that our best fits are still com- 
patible with the already mentioned upper limits on the 
photon f ract i on o f the observed ultrahigh energy cosmic 
We observe that the value of the neu- 



rays U2, 155, n5J 



trino mass found in both the halo background scenario, 



with ; 



3.71 



+ 1.40(3.27) 



eV, as well as in the extragalac- 



77 +0-48(1.36) y 
u " 1 1 -0.30(0.51) c v ' 



-1.12(1.96) 

tic background scenario, with m v 
are compatible, at about the 2 a level, with the corre- 
sponding values found in the case of strong UHE7 atten- 
uation. As before, in the "no" background scenario, we 
find a strong dependence of the fitted mass on the value 

Of -Emin- 

Besides the no UHECR background scenario, the Z- 
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FIG. 8: The available UHECR data with their error bars and 
the Xmin = t4 fits (solid lines) from Z-bursts, for a halo (top), an 
extragalactic (middle), and no (bottom) UHECR background 
(long-dashed lines) , and a strong UHE7 attenuation such that 
the Z-burst photons can be neglected and only the Z-burst 
protons (dash-dotted lines) have to be taken into account in 
the fit (a = 0, h = 0.71, fl M = 0.3, fi A = 0.7, z max = 2). 
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FIG. 9: The available UHECR data with their error bars 
and the best fits ffomZ-bursts, for an energy attenuation of 
photons as in Ref. [116 and a vanishing EGMF (a = 0, h = 
0.71, Q M = 0.3, Q.a = 0.7, z max = 2). Top: Best fit for the 
case of a halo UHECR background (solid line), corresponding 
to the sum of the background protons (long-dashed), the Z- 
burst protons (dash-dotted) and the Z-burst photons (short- 
dashed) . Middle: The case of an extragalactic UHECR back- 
ground. Bottom: The case of no UHECR background above 
log(£ min /eV) = 19.7. 



TAB LE I I: Results of fits, for a photon attenuation as in 
Ref. Jll§ ("minimal" URB) and vanishing EGMF (h = 

2). Top: Assuming a 
"1). Middle: 



0.71, = 0.3, Q A = 0.7,z m ax = 2). Top: 

halo UHECR background according to Eq. (| 

Assuming an extragalactic UHECR background according to 
Eq. (B5J). Bottom: Assuming no UHECR background above 
£min, for different values of the lower end .Emm of the fit 
(q = 0). 



Halo UHECR background -I 


- minimal URB 




a 


m„ [eV] 


Xmin 


A 


B 


P 


-3 

3 


4 46+1-78(3.99) 
^•^ u -1.36(2.39) 
o ?1 +1.40(3.27) 
-1.12(1.96) 
o m +1.13(2.62) 
,5 - Ui -0.92(1.56) 


15.69 
15.36 
15.20 


4.76 ■ 10 43 
1.05 ■ 10 44 
1.55 ■ 10 44 


1060 
827 
642 


3.146 
3.166 
3.171 


EG UHECR background + minimal URB 


ct 


m„ [eV] 


Xmin 


A 


B 




-3 

3 


n YQ +0.49(1.40) 
u ' 1 a -0.31(0.53) 
„ 77 +0.48(1.36) 
U ' ' ' -0.30(0.51) 
„ „ q +0.45(1.32) 
u -' a -0.31(0.53) 


24.19 
24.52 
24.77 


5.73 ■ 10 31 
6.44 ■ 10 31 
7.55 ■ 10 31 


125 
120 
122 


2.466 
2.467 
2.468 


No UHECR background + 


minimal URB 




log(EUi»/eV) 


m„ [eV] 


Xmin 


A 


B 


P 


19.4 
19.5 
19.6 
19.7 
19.8 
19.9 
20.0 


„ n „+0.78(1.68) 
J - UO -0.58(1.05) 

ce+0.76(1.68) 
Z ' JO -0.55(0.97) 
9 ,-,+0.84(1.87) 
Z - J1 -0.60(1.05) 

1 n, +0.75(1.76) 
1 -° i -0.51(0.86) 
9 n .+0.95(2.17) 
z ' uy: -0.67(l.ll) 
o n , +1.34(3.03) 

^•^-l.oea^s) 

. „ 7 +2. 12(4.60) 
^■ ut -1.74(2.95) 


14.43 
11.88 
11.86 
8.02 
7.74 
4.93 
2.33 




1380 
1111 
1079 
680 
812 
1699 
4881 





burst determination of the neutrino mass seems reason- 
ably robust. Its robustness with regard to changes of 
presently unknown quantities within t heir anticipated 



III 



and [V 



and 



variations is further illustrated in Tables 
in Fig. For a wide range of cosmological source evo- 
lution (a = —3 -j- 3), Hubble parameters h = 0.61 -r 0.9, 
Qm , ^A, z max = 2-^5, for variations of the possible relic 
neutrino overdensity in our GZK zone within the limits 
discussed in Section II D, and for different assumptions 



about the diffuse extragalactic photon background, the 
results remain within the above error bars. The main un- 
certainties concerning the central values originate from 
the different assumptions about the background of ordi- 
nary cosmic rays. In the case that the ordinary cosmic 
rays above 10 18 5 eV are protons and originate from a 
region within the GZK zone of about 50 Mpc ("halo"), 
the required mass of the heaviest neutrino seems to lie 
between 2.1 eV < m v < 6.7 eV at the 68 % C.L. (a < 0, 
cf. below), if we take into account the variations between 
the minimal and moderate URB cases and the strong 
UHE7 attenuation case (cf. Fig. [To]). Note, that a value 
of m y — 0.07 eV, as studied recently in Ref. |56| in a sim- 
ulation of the Z-burst scenario in the context of a halo 
background model, seems to lie more than 3 a away from 
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TABLE III: Results of fits, for a photon attenuation as in 
Fig. |§] (bottom), short-dashed line ("moderate" URB) and 
vanishing EGMF (h = 0.71, Sl M = 0.3, Q A = 0.7,z max = 2). 
Top: Assuming a halo UHECR background according to 
Eq. ([mJ). Middle: Assu ming an extragalactic UHECR back- 
ground according to Eq. (Ba). Bottom: Assuming no UHECR 
background above i?mm, for different values of the lower end 
Em in of the fit (a = 0). 



Halo UHECR background + moderate URB 


a 


m v [eV] 


Xmin 


A 


B 


P 


-3 

3 


a or + l. 97(4.31) 
45(2. 57) 
q fi9 +1.61(3.66) 
°-° -1.20(2.13) 
9 04+1-30(3.00) 
-1.01(1.72) 


15.69 
15.44 
15.26 


2.16 ■ 10 43 
2.38 ■ 10 43 
4.91 ■ 10 43 


848 
662 
514 


3.112 
3.113 
3.149 


EG UHECR background + moderate URB 


a 


m v [eV] 


Xmin 


A 


B 


P 


-3 

3 


n q 9 +0.29(0.84) 
U-J^-0. 17(0.27) 
„ q 9 +0.28(0.80) 
U - 3Z -0.17(0.26) 
„ q, +0.28(0.78) 
U " 31 -0.17(0.26) 


21.15 
21.29 
21.37 


4.95 ■ 10 31 
5.40 ■ 10 31 
5.94 • 10 31 


60 

58 
56 


2.465 
2.466 
2.466 


No UHECR background + 


moderate URB 




log(£U n /eV) 


m v [eV] 


Xmin 


A 


B 


P 


19.4 
19.5 
19.6 
19.7 
19.8 
19.9 
20.0 


9 „i +0.74(1.65) 
Z -° -0.59(1.05) 
„ , „+0.73(1.65) 
Z ' iO -0.52(0.92) 

9 n/i +0 - 81(1 - 84) 

Z - U -0.56(0.97) 
1 , 9 +0.59(1.50) 
i ' iZ -0.40(0.64) 

1 nzi +0 - 75(1 - 92) 

1 - u, -0.44(0.69) 
i RR +1.31(3.03) 
1 - OO -0.81(1.23) 
„ oi+2.15(5.03) 
Z - 0i -1.44(2.24) 


15.01 
11.99 
11.82 
5.53 
5.49 
4.26 
2.75 




1087 

820 

755 

369 

341 

652 

1728 





our best fit values. The much more plausible assumption 
that the ordinary cosmic rays above 10 18 5 are protons of 
extragalactic origin leads to a required neutrino mass of 
0.08 eV < m v < 1.3 eV at the 68% C.L. {a < 0) (cf. 
Fig. 0). 

We performed a Monte Carlo analysis studying higher 
statistics. I n th e near future, the Pierre Auger Obser- 



vatory 161 , 162] will provide a ten times higher statis- 
tics, which reduces the error bars in the neutrino mass 
to about one third of their present values. 

Let us consider next in more detail the 7 ray spectra 
from Z-bursts, notably in the ~ 100 GeV region. As 
illustrated in Fig. [Uj, the EGRET measurements of the 
diffuse 7 back grou nd in the energy range between 30 MeV 
and 100 GeV 163 put non-trivial constraints on the evo- 
lution parameter a. Whereas different values of a lead 
to very similar spectra in the ultrahigh energy region 
- which determines the neutrino mass - they are easily 
discriminated in < TeV photons. Only evolution param- 
eters a < seem to be compatible with the EGRET mea- 
surements (cf. Fig. [ll]) , quite independently of different 
assumptions about the URB. These numerical findings 
are in fairly good agreement with other recent simula- 
tions fjl. 



TABLE IV: Results of fits, for strong and minimal photon at- 
tenuation (cf. Tables j] and |n|, respectively) , for extremal val- 
ues of the (reduced) Hubble constant (a = 0, Qm = 0.3, Q,a = 



0.7, 2„ 



= 2). 



EG UHECR background + strong UHE7 attenuation 


h 


m„ [eV] 


Xmin 


A 


B 


P 


0.61 
0.90 


„ 1 q +0.19(0.62) 
-0.12(0.17) 
„ 91 +0.20(0.64) 
U ' Z -0.11(0.18) 


28.27 
24.20 


1.47 ■ 10 31 
8.93 • 10 32 


126 
168 


2.452 
2.529 


EG UHECR background + minimal URB 


h 


m„ [eV] 


Xmin 


A 


B 


P 


0.61 
0.90 


n fi4+ () - 42 (i' 2 3) 

u,u -0.25(0.44) 
1 n9 +0.64(1.70) 
1 - uz _0.39(0.65) 


26.28 
22.39 


1.46 • 10 31 
1.20 ■ 10 33 


97 
162 


2.452 
2.531 


EG UHECR background + moderate URB 


h 


m v [eV] 


Xmin 


A 


B 


P 


0.61 
0.90 


„ 9fi +0.25(0.71) 
u ' zu -0.14(0.22) 
n 4n +0.37(1.03) 
u -0.20(0.32) 


22.91 
19.53 


8.03 • 10 30 
8.41 • 10 32 


52 
64 


2.435 
2.529 



The necessary UHEC^ flux at the resonant energy E]f s 
is obtained via Eq. (p37j) from the fitted overall normal- 
ization B, 



CT-os pros = g . 1Q - 



1 



D 



m 2 s sr R„ 



Hi 



S-hB, 



(40) 



where Dh = c/Hq = 3 • 10 3 Yr 1 Vfpc is the present size 
of the universe and where, on account of Eqs. (|2^) and 
( p8| ) , the ratio of the maximally considered distance of Z 
production, i? max , and the size of the universe is given 
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FIG. 10: Summary of the masses of the heaviest neutrino 
required in the Z-burst scenario, with their 1 a (solid) and 2 a 
(dotted) error bars, for the case of an extragalactic and a halo 
background of ordinary cosmic rays and for various assump- 
tions about the diffuse extragalactic photon background in the 
radio band (a = 0, h = 0.71, fi M = 0.3, tt A = 0.7, z max = 2). 
From left: strong 7 attenuation, moderate and minimal URB. 
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FIG. 11: The available UHECR data with their error 
bars and the best fit from Z-bursts, for various cosmolog- 
ical evolution parameters a and an energy attenuation of 
photons as in Fig. |^ (bottom) exploiting a "minimal" URB 
(h = 0.71, Q,m ~ 0.3, Q,a = 0.7, z max = 2). Also shown is the 
diffuse 7 background in the energy ra nge between 30 MeV 
and 100 GeV as measured by EGRET pi (solid). 



by 



Dh 



dz ■ 



(i + z) y/n M (i + z f + n A 



(41) 



This ratio is of order one; it equals 0.73 for Qm = 0.3, 
SIa = 0.7 and z max = 2. The required fluxes are sum- 
marized in Fig. 
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together with some existing upper 
limits and projected sensitivities of present, near future 
and future observational projects. They appear to be 
well below present upper limits and are within the ex- 



pected sensitivity of AMANDA |L64], [16§, RICE p8[ , 
and Auger 166]. Clearly, these fluxes are higher than the 
ones advocated in Ref. [|52f based on local neutrino over- 
densities /„ = 10 2 -j- 10 4 on scales of about 5 Mpc. How- 
ever, since we also took into account a background from 
ordinary cosmic rays, our normalization of the Z-burst 
component is different and correspondingly our fluxes are 
somewhat less than a factor of /„ higher. They are con- 
sistent with the ones found in Refs. [^5|, [36|. As far as 
power-law extrapolations of the fluxes, F v oc E" 1 , below 
the resonance energy are concerned, we find, in agree- 
ment with Refs. |5^, [33| ^35|, t hat i ndices 7>1.5 are ex- 
cluded by the Fly's Eye limits fL2§ (cf. Fig. 0). 

The required neutrino flux, for power-law extrapola- 
tions with index 7>1 to lower energies, is larger than 
the theoretical upper limit from "hidden" hadronic as- 
trophysical sources, from which only photons and neutri- 
nos can escape (cf. Fig. |il("MPR"j). From an analysis of 
the assumptions behind this limit [127| one finds that one 
has to invoke hidden astrophysical sources whose photons 
somehow don't show up in the diffuse 7 ray background 
measured by EGRET |i4], [35|. In general, astrophys- 
ical sources for the required UHE neutrinos should be 



FIG. 12: Neutrino fluxes, F = | Y%=x( F »i + F ^)> recmired 
by the Z-burst hypothesis for the case of a halo and an ex- 
tragalactic background of ordinary cosmic rays, respectively 
(a = 0, h = 0.71, Q, M = 0.3, Q. A = 0.7, z max = 2). Shown are 
the necessary fluxes obtained from the fit results of Table |] 
for the case of a strong UHE7 attenuation. The horizontal 
errors indicate the 1 a (solid) and 2 a (dotted) uncertainties 
of the mass determination and the vertical errors include also 
the uncertainty of the Hubbl e exp ansion rate. Also shown are 
upper limits from Fly's Eye |123[ on _F„ C + Fa e and the Gold- 



stone lunar ultrahigh energy neutrino experiment GLUE [ 124 
on X/ Q =e S Fv " "t ^fa)> as weu as projected sensitivities of 



166 1 



AMANDA (16J, |165| on F V)1 + F^ and Auger 
F„ e + Fu. . Th e sensitiviy of RICE is comparable to the one 
of Auger Egl. 



distributed with a < 0, accelerate protons up to en- 
ergies > 10 23 eV, be opaque to primary nucleons and 
emit secondary photons only in the sub-MeV region. It 
is an interesting question whether such challenging con- 
ditions can be realized in BL Lac objects, a class of 
active galactic nuclei for which some evidence of zero 
or n egati ve cosmological evolution has been found (see 
Ref. [167 1 and references therein) and which were recently 
discussed as possible sources of the highest energy cos- 
[|168|. Alternatively, one may invoke top-down 
for the sources of the highest energy cos- 
mic neutrinos such as u nstable superheavy reli c particle 
decays ||| H Ell EM EH ' EH Ell EH or to P " 




logical defect decays |147[ 

It should be stressed that, besides the neutrino mass, 
the UHECi/ flux at the resonance energy is one of the 
most robust predictions of the Z-burst scenario which 
can be verified or falsified in the near future. In this 
connection it is worthwhile to recall (cf. Sect ion IIP ) 



that the current limits on neutrino degeneracies [153, 154 
do not allow a remarkable uniform enhancement due to 
lepton asymmetries which otherwise would be welcome 
for a relaxation of the huge flux requirement pm . 
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IV. DISCUSSION AND CONCLUSIONS 

We have presented a detailed comparison of the pre- 
dicted spectra from Z-bursts - resulting from the reso- 
nant annihilation of ultrahigh energy cosmic neutrinos 
with relic neutrinos - with the observed ultrahigh en ergy 
cosmic ray spectrum, extending our earlier study |pl| . 
The mass of the heaviest relic neutrino turned out to 
have to lie in the range 2.1 eV < m„ < 6.7 eV 
on the 1 a level, if the background of ordinary cosmic 
rays above 10 185 eV consists of protons and originates 
from a region within a distance of about 50 Mpc. In 
the phenomenologically most plausible case that the or- 
dinary cosmic rays above 10 18 ' 5 are protons of extra- 
galactic origin one is lead to a required neutrino mass 
in the range 0.08 eV < m v < 1.3 eV at the 68% C.L.. 
We have also analysed a case where there is no back- 
ground of ordinary cosmic rays above, say, 10 19 ' 7 eV. 
Here, we find 0.24 eV < m v < 2.6 eV. The above 
neutrino mass ranges include variations in presently un- 
known quantities, like the amount of neutrino cluster- 
ing, the universal radio background, and the extragalac- 
tic magnetic field, within their anticipated uncertainties. 
If it turns out that the highest energy cosmic rays are 
extragalactic protons rather than photons one has to in- 
voke a large universal radio background and/or large ex- 
tragalactic magnetic field to suppress the photons from 
Z-bursts. In this case and for our extragalactic back- 
ground scenario, the neutrino mass range narrows down 
to 0.08 eV < m„ < 0.40 eV, with a best fit value of 
m v = 0.20 eV. 

It is remarkable, that the mass ranges required in 
the Z-burst scenario coincide nearly perfectly with the 
present knowledge about the mass of the heaviest neu- 
trino from oscillations and tritium j3 decay from Eqs. (|^) 
and (||), 0.04 eV < m„ 3 < 2.5 eV, in a three flavour or, 
from Eqs. (@) and (§), 0.4 eV < m Vi < 3.8 eV, in a four 
flavour scenario. Our values, in the extragalactic back- 
ground scenario, are still compatible with a hierarchi- 
cal neutrino mass spectrum, with the largest mass sug- 
gested by the atmospheric neutrino oscillation, m U3 ss 



with a projected sensitivity of mp > 0.3 eV within 6 



0.04-^0.08 eV. However, our best fit values 



point more into the direction of a quasi-degenerate sce- 
nario, m Ui rj rriy ^> \J Am 2 tm . It is amusing to observe 
that the recently reported evidence for neutrinoless dou- 
ble beta decay, with 0.11 eV < (m„) < 0.56 eV, if true, 
would point also into the same direction |Q . 

The above neutrino masses are in a range which 
can be explored by near- future laboratory experi ments , 
like the KATRIN tritium j3 decay experiment [169|, 



7 years, or the Oz/2 /3 decay experiments NEMO-3 [17C] 
and CUORE fll7l]] , which aim at (m u )>0.l eV, and 
their next- gene ration follow-up projects GENIUS [172] 
and EXO [173|, with (m„)>0.01 eV. Our mass range 
also compares favourably with the expected sensitiv- 
ity Yli(9vt/2) m vi > 0.3 eV of near-future global cos- 
mological analyses involving also new CMB measure- 
ments |H [174] , as well as mp > 0.75 4- 1.8 eV ex- 
pected from neutrino observations from supernovae ex- 



plosions Il75[ 176 



From our fits, we could also determine the necessary ul- 
trahigh energy neutrino flux at the resonance energy. Its 
prediction, in the context of the Z-burst scenario, turned 
out to be of similar robustness as the prediction of the 
mass of the heaviest neutrino. It was found to be con- 
sistent with present upper limits and detectable in the 
near future by the already operating neutrino telescopes 
AMANDA and RICE, and by the Pierre Auger Obser- 
vatory, presently under construction. A search at these 
facilities is the most promising and timely step in test- 
ing the Z-burst hypothesis. One does not have to wait 
for future projects such as EUSO [177| for this critical 
investigation. 

The required neutrino fluxes are enormous. If such 
tremendous fluxes of ultrahigh energy neutrinos are in- 
deed found, one has to deal with the challenge to explain 
their origin. It is fair to say, that at the moment no 
convincing astrophysical sources are known which meet 
the requirements for the Z-burst hypothesis, i.e. which 
have no or a negative cosmological evolution, accelerate 
protons at least up to 10 23 eV, are opaque to primary nu- 
cleons and emit secondary photons only in the sub-MeV 
region. Ideas to create the required neutrino fluxes from 
the decays of superheavy particles or topological defects 
do not seem to be too appealing since they invoke further 
new physics beyond neutrino masses. 
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